Dynamical instability and domain formation in a spin-1 Bose condensate 
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We interpret the recently observed spatial domain formation in spin-1 atomic condensates as a 
result of dynamical instability. Within the mean field theory, a homogeneous condensate is dy- 
namically unstable (stable) for ferromagnetic (antiferromagnetic) atomic interactions. We find this 
dynamical instability naturally leads to spontaneous domain formation as observed in several recent 
experiments for condensates with rather small numbers of atoms. For trapped condensates, our 
numerical simulations compare quantitatively to the experimental results, thus largely confirming 
the physical insight from our analysis of the homogeneous case. 
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Spatial domains or pattern formation is a common fea- 
ture of nonlinear dynamics in extended systems. It has 
been actively researched in nonlinear optics 0, classical 
fluids 0, granular materials 0, and recently in atomic 
Bose-Einstein condensates (BECs) 0, S S 0, B H E3 • 
ft is generally understood that the unstable modes of 
a dynamically unstable system can grow exponentially 
and eventually lead to the appearance of spatial domain 
structures that last for a long time. 

Many earlier studies have suggested interesting mech- 
anisms for spontaneous domain formation in atomic con- 
densates 0, IM IE • Most focus on the single- or two- 
component condensates, where the number(s) of atoms 
for each component is conserved. The dynamical insta- 
bility due to attractive atomic interactions is the most 
prominent among all proposed scenarios for domain for- 
mation 0, IE 13 ■ The attractive interaction in a single- 
component condensate is also believed to be responsible 
for the formation of a train of solitons, consistent with 
the fact that it is dynamically unstable Q. For a two- 
component condensate, again it is found that an effective 
attractive interaction is responsible for the dynamical in- 
stability and domain formation 0, @, 0] ■ 

Several groups have also studied three-component, or 
spin-1 condensates (F = IV which are distinct as the spin 
mixing interaction [q. llCtlll| allows for exchanging atoms 
among spin components 2\mp = 0) <-> \mp = +1) + 
\nip = —1) (hereafter as |0), |+), and |— )). The number 
of atoms for each component therefore can change, but 
the total number of atoms and the system magnetization 
are conserved. Significant interest now exists for spin-1 
condensates because of the recent progress from several 
experimental groups |12|. in particular, the observation 
of spontaneous domain formation in 87 Rb condensates 
|13|. Robins et al. were among the first to study dy- 
namical instability in a spin-1 condensate [ljj. They dis- 
covered a particular type of stationary state dynamically 
unstable for ferromagnetic interactions, evidenced by the 
sudden collapse when propagated with Gross-Pitaevskii 
(GP) equations, presumably resulting from the amplifi- 
cation of numerical discretization errors. Through exten- 
sive numerical simulations, Saito and Ueda also investi- 
gated very recently the spontaneous multi-domain for- 



mation induced by the dynamical instabilit y in a spin-1 
condensate with ferromagnetic interactions |ll| . A clear 
picture, however, is still lacking as indicated by the gen- 
eral lack of comparisons with experimental reports. Our 
work aims at providing a complete understanding for do- 
main formation in a spin-1 condensate. 

To begin with, we consider a homogeneous condensate 
at an off-equilibrium state initially. For example, a spin-1 
condensate in the ground state at a certain nonzero mag- 
netic (B-) field for time t < will become off-equilibrium 
when the external B-field is changed for t > 0. This 
causes the spin-1 condensate to collectively oscillate anal- 
ogous to a nonrigid pendulum as we recently showed . 
In addition, we assume the condensate size is much larger 
than the spin healing length at least in one direction so 
that domains may be formed. Within our mean field 
description, the evolution of a spin-1 condensate is de- 
scribed by the coupled GP equations [15| 
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where Ti = — (ft 2 /2m)V 2 + V cyiii + c n, <&j is the j-th spin 
component condensate wave function, and rij = |$j| . 
c = 4nh 2 (a + 2a2)/3m and c 2 = 47r?i 2 (a 2 — ao)/3m 
with ao and <z 2 the scattering lengths for the two collid- 
ing atoms in the symmetric channels of total spin and 
2, respectively. The interaction is ferromagnetic (antifer- 
romagnetic) if c 2 < (> 0). 
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due to the conservation of atomic density (n = n + + uq + 
n_) and the magnetization (m = n + — Equations (J2J) 
defines an energy conserving dynamics, with the effective 
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energy per unit of volume given by 
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We note Eq. for a homogeneous condensate differs 
from a trapped one even under single spatial mode ap- 
proximation despite sharing the same dynamical Eq. J2J. 

Within the mean field approximation, the average spin 
of a condensate / = f x x + f y y + mz, where fj = (Fj) 
with F X , V:Z being the spin-1 matrices, is also conserved in 
addition to the conservations of n and m [l6l | . The energy 
functional Eq. © thus becomes £ = \c§n 2 + hc 2 f 2 , if 



We note that the mean field theory 
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model cannot be applied to extreme cases such as A^o 
and No = N, where quantum effects are important. 

We adopt three approaches to study dynamical stabil- 
ity of the off-equilibrium collective oscillations of a con- 
densate: the effective potential method; the Bogoliubov 
method; and direct numerical simulations. By going into 
a rotating frame, an entire orbit reduces to a station- 
ary point in the phase space 01 ■ The effective potential 
then becomes T = (c /2)n 2 + (c 2 /2) (m 2 + / 2 + / 2 ) - 
fin — nm — d x f x — S y f y , where parameters n, 6 X , 5 y } = 
{con, c 2 m, C2/Z, C2/2,} defines the rotating frame and are 
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Our system is dynamically stable if its Hessian matrix 
of T with respect to {n, m, f x , f y } is positive definite and 
dynamically unstable if the Hessian matrix has any neg- 
ative eigenvalue. It is easy to check that the eigenvalues 
of the Hessian matrix are {co, 02,02,02}. Thus an antifer- 
romagnetically interacting spin-1 condensate is dynami- 
cally stable, while a ferromagnetically interacting one is 
dynamically unstable since C2 < 0. 

We next employ the Bogoliubov transformation to find 
out the corresponding unstable modes. Starting from the 
stationary point in the rotating frame as found above, the 
equation of motion for collective excitations can be cast 
in a matrix form |l7^ as A4 ■ x — Kujx, with a vector 
x = (S^+, 5^o, d^-,5^l,5^,5^*_) T . SVj and SV* de- 
note the deviations from the stationary point, and the 
associated matrix is 
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e k = fr 2 k 2 /2m is kinetic energy of the collective excitation 
mode with wave vector k. 

The eigen-frequencies of the Bogoliubov excitations 
are obtained from the characteristic equation det(.M — 
HujI) = 0, explicitly given by 
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The corresponding modes are termed as density modes 
(solid lines), spin modes (dotted lines), and quadrupolar 
spin modes (dashed lines as in Fig. GJ by Ho [13. Fi gurc 
H^a) and Gib) show the real and imaginary parts of the 
typical dispersion relation for a 87 Rb spin-1 condensate 
|18| . respectively. All frequencies are real for antiferro- 
magnetically interacting (e.g. 23 Na) condensate. 

Our analysis here parallels that of Refs. 0, || for 
a two-component condensate. We find two interest- 



FIG. 1: (Color online) Real part (a) and imaginary part (b) 
of a typical Bogoliubov spectrum for a spin-1 condensate 
with ferromagnetic interaction. Panel (c) shows the small- 
est wavelength of a homogeneous spin-1 87 Rb condensate at 
n = 1.9 x 10 14 cm -3 . 



ing features in Fig. ^b). One of them at fc_ is the 
most unstable mode at the maximum imaginary fre- 
quency, it determines short time behavior such as the 
time scale for domains to emerge. The other is the 
largest wave vector k m with an (infinitesimal) imagi- 
nary frequency which determines the long time behav- 
ior such as the final domain size. From Eq. JSJ, we 
find the time scale for the emergence of a domain is 
~ l/|tj(fc_)| = hj (|c s | Vl — f 2 ) and the domain width 



is about A TO = 2n/k m = h/y2m\c s \(l + \/l — f 2 ). It is 
typically of the order of spin healing length. Figure ^c) 
displays the f-dependence of A m at n — 1.9 x 10 14 cm -3 
for a homogeneous 87 Rb spin-1 condensate. For the do- 
mains to form, a condensate has to be larger than A m , at 
least in one direction. 

We note the spin domain formation as discussed here 
is different from striation patterns as observed in (an- 
tiferromagnetic) 23 Na condensates. The stripe patterns 
arise from interplay of an external B-field, a field gradi- 
ent, and immiscibility among different spin components 
No domains were observed in Stenger et aL's ex- 
periment at negligible B-fields where the |+) and |— ) 
components coexist. At finite values of B-fields, phase 
separation between the |0) and the |+/— ) components 
occurs [l9j . In Miesner et al. and Stamper-Kurn et a/.'s 
experiments 6] , only two spin components were involved 
due to the relatively large bias B-field (~ 15 G). 

Our analysis shows the formation of spin domains is a 
direct consequence of dynamic instability for a conden- 
sate with ferromagnetic interactions. To provide a clearer 
physical picture for domain formation, we now work in 
the lab frame. We focus on d£/dm, which in fact cali- 
brates the formation of the spin domain. We find 
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FIG. 2: (Color online) Surfaces of dS/dm = (top) and 
the cross-section at 9 — (bottom). The plus signs denote 
d£/dm > and the minus signs denote d£/dm < 0. 



derivative is zero. The region below the saddle surface in 
Fig. [SJof an orbit is unstable if c-i < 0. Here the meaning 
of "unstable" is generalized, referring to the dynamical 
property where the local magnetization tends to deviate 
further from m — 0. For example, in the lower right al- 
lowed region of d£/dm < 0, Am > is required to lower 
the local energy £ w £ (m) + (d£/dm)Am. Thus m tends 
to increase. Similarly the lower left allowed region would 
make m decrease. The combined effect is dynamically 
unstable orbits, separation of |+) and |— ) component, 
and the eventual formation of spin domains. 

For antiferromagnetic interactions (c~2 > 0), 8 usually 
oscillates around n. Thus d£/dm > for m > and 
d£ I dm < for m < 0. So the magnetization always os- 
cillates around zero and no domain forms. This coincides 
with the findings of the dynamical stability analysis for 
an antiferromagnetically interacting condensate. 

Finally we perform numerical simulations of Eq. Q to 
confirm the mechanism of dynamical instability-induced 
spin domain formation. The initial conditions are as in 
the experiment ^sjj with 87 Rb condensates (Nq(0)/N — 
0.744, 6»(0) = 0, for the ground state of N = 2.0 x 10 5 
at B — 0.3 G and M = 0), in a trap Kxt(^) = 
(m/2)(w 2 2! 2 + uo 2 y y 2 + uj 2 z z 2 ) with uj x = uj y = (2tt)240 Hz 
and co z — (27r)24 Hz. In one of the simulations, we inten- 
tionally include additive small white noise (~ 1.0 x 10 -5 ), 
although still much larger than numerical errors |20j dur- 
ing the propagation. We find that it takes a shorter time 
for the |+) and |— ) components to separate when white- 
noise is included. Figure |3] shows the evolution of ax- 
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FIG. 3: (Color online) Typical evolutions for spin domain for- 
mation in a 87 Rb condensate. The initial state is the ground 
state at B = 0.3 G. The B-field is then set to zero and a 
small white noises are added throughout the evolution. The 
left contour plot is for the |+) component. The right col- 
umn shows the density distribution of all three components at 
times t — 0,160,320 (1/oj z ). Solid, dash-dotted, and dashed 
lines denote respectively the |+), |0), and |— ) component. 
Dotted lines are for the total density. The axial density is 
n z = J \ 2 2-Krdr. a z — ^/h/mu> z ~ 2.2/^m and the 

average condensate density is ~ 1.9 x 10 14 cm -3 . 



ial density distributions. Phase separation between the 
|+) and |— ) components is seen accompanied by the for- 
mation of domains. This proves again that dynamical 
instability causes the formation of domains. The do- 
main width (an upper limit) as estimated from X m = 

h/ ^J2m\c2\(n)(l + \J\ — W) w 15/im, is consistent with 
both simulations and experimental observations |13|. 

In a spinor condensate, spin wave excitations normally 
refer to dynamically stable (or relatively more stable) col- 
lective modes. Once excited, they lead to coherent cyclic 
dynamics in both spatial and temporal dimensions. Spin 
domains, on the other hand, refers to unstable modes, 
with a fixed pattern in the long time limit. 

In conclusion we have presented a systematic study 
of dynamical stability and the accompanied mechanism 
for domain formation in a spin-1 condensate. Our re- 
sults affirm that a ferromagnetically interacting conden- 
sate is dynamically unstable and evolves spontaneously 
into multi-domain structures contrary to dynamically 
stable antiferromagnetic condensates. Our work provides 
a clear physical picture for recently observed spontaneous 
domain formations in spin-1 condensates [l3|. 

This work is supported by NASA and NSF. 
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